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ABSTRACT 



The Telescope Array (TA) collaboration has measured the energy spectrum 
of ultra-high energy cosmic rays with primary energies above 1.6 x 10 18 eV. This 
measurement is based upon four years of observation by the surface detector 
component of TA. The spectrum shows a dip at an energy of 4.6 x 10 18 eV 
and a steepening at 5.4 x 10 19 eV which is consistent with the expectation from 
the GZK cutoff. We present the results of a technique, new to the analysis of 
ultra-high energy cosmic ray surface detector data, that involves generating a 
complete simulation of ultra-high energy cosmic rays striking the TA surface 
detector. The procedure starts with shower simulations using the CORSIKA 
Monte Carlo program where we have solved the problems caused by use of the 
"thinning" approximation. This simulation method allows us to make an accurate 
calculation of the acceptance of the detector for the energies concerned. 



1. Introduction 



One of the most powerful tools for studying the origin of ultra-high energy cosmic rays 
(UHECRs) is their energy spectrum, which manifests several features that reveal important 
information about the cosmic rays, their sources, and their propagation across cosmological 
distances. One example is the h i gh-ene rgy (4 — 6 x 10 19 eV) suppression in the spectrum 



which was predicted by 



Zatsepin &: Kuzmin 



(Il966f ). and is called 



Greisenl (119661 ) and by 

the GZK cutoff. These authors predicted a strong suppression in the spectrum due to the 
interaction of cosmic rays with photons of the cosmic microwave background radiation. 
They pointed out that a spectrum suppression is expected in both cosmic protons (by 
photo-pion production) and heavier nuclei (by spallation) for particles traveling more than 
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50 Mpc from their sources. If cosmic rays are protons there should also be a dip in the 
spectrum, ca used by e + e~ pair production i n the same interactions, at an energy of about 



5 x 10 eV (iBerezinsky fc Grigor'eva 



19881 ) . For heavier nuclei, interactions with the 



background photon flux do not cause such a dip. 



The AGASA experiment (ITakeda et al.lll998l . 120031 ). comprised of a surface array of 
111 scintillation counters, was the first detector to be large enough to test this theory with 
sufficient statistics. However, they did not observe the suppression. The first experiment 
to observe the GZK cutoff was the High Resolution Fly's Eye (HiRes) experiment 



(lAbbasi et al. 



20081 ). which consisted of fluorescence detectors located atop two desert 



mountains in western Utah. HiRes reported a cutoff energy of (5.6 ± 0.5 ± 0.9) x 10 19 eV, 
which is consistent with a suppression of protons. They also observed the ankle structure: 
a hardening of the spectrum, at an energy of 4.5 ± 0.05 ± 0.8 x 10 18 eV as expected for 
cosmic protons. (For both of these values, the first uncertainty is statistical and the second 
systematic.) HiRes also published measurements of the shower maximum s lant depth 
(I mm ) that indicated a predominately light composition above 2 x 10 18 eV (lAbbasi et al. 



2010|). 



A somewhat different picture is seen by the Pierre Auger Observatory (PAO), located in 
Argentina. The PAO consists of a surface detector (SD) of 1600 water tanks, accompanied 
by 24 fluorescence telescopes which are equally apportioned between four sites located at the 
SD corners. The PAO also obser ves the high-energy suppression, but at (2.9 ± 0.2) x 10 19 



eV ( [Abraham et al. 



2008 



2010bl ). They see the ankle also, b ut their X mnx results m ay 



2010ah . One 



indicate that the composition is heavy at the highest energies (1 Abraham et al.l 
possible interpretation of the PAO results is that the high-energy suppression is caused by 
spallation of heavy nuclei. The cause of the ankle would need to be explained by a separate 
mechanism. 
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The Telescope Array (TA) experiment, also in western Utah, is the largest experiment 
studying ultra-high energy cosmic rays in the northern hemisphere. A layout of the TA 
experime nt is shown in Figure U It consists of a surface detector of 507 scintillation 



counters (lAbu-Zayvad et al. 



Abu-Zayyad et al. 



2012 



2012q ). plus 38 fluorescence telescopes (ITokuno et al. 
2012bl ) located at three sites overlooking the SD. 24 of the telescopes, 
deployed at the two southern sites, were built new for the experiment, and the 14 telescopes 
at the northern site are reconditioned HiRes telescopes. TA combines the experimental 
techniques of AGASA and HiRes, in order to understand the difference between their 
results. 

This paper reports on a measurement of the cosmic ray spectrum above 1.6 x 10 18 eV 
made by the TA SD over approximately four years of observation between May 11, 2008 
and May 20, 2012. For this study, we used an analysis method that while standard for 
fluorescence detectors, is being successfully implemented for the first time for a surface 
array studying cosmic rays in the ultra-high (> 1 EeV) energy regime. Instead of restricting 
our analysis to a domain where we expect 100% efficiency, as has been done by previous 
surface detector experiments at these energies, the TA SD detector aperture is calculated 
using extensive a ir showers generated in detail by the CORSIKA simulation package 



Heck et al. 



19981 ). accompanied by a full GEANT simulation of the detector (lAllison et al. 
20061 ). Another important aspect of this technique, new to surface detectors operating in 
the ultra-high energy regime, is the validation of the simulation by comparisons of key 
distributions from the data to those obtained from the Monte Carlo (MC) simulation. 
Moreover, our study overc omes the inab ility of "thinned" simulated showers (e.g. as used 



in CORSIKA and AIRES flSciutto 



2002])) to reproduce the particle density and arrival time 



fluctuations far from the core. The solution app 



replicates a non-thinned simulation (jStokes et al 
most of the detector data is collected. 



ied is a novel dethinning technique that 



20121 ) at the lateral core distances where 
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The TA Surface Detector 



Each counter of the TA SD consists of two layers of 1.2 cm thick plastic scintillator, 
both 3 m 2 in area. Photons produced by ionizing particles passing through the counters 
are collected by wavelength shifting fibers and read out by photomultiplier tubes, one for 
each layer. A histogram of pulse heights, triggered by a coincidence between the two layers 
within an individual SD, is collected every 10 minutes. This histogram is dominated by 
single muons with a count rate of ~ 700 Hz. Each 10-minute histogram is used to calibrate 
the associated scintillator to the pulse height of a minimum ionizing particle (MIP) with 
~ 1% accuracy. The SD array trigger requires at least three adjacent counters with pulse 
heights over 3 MIP to fire within 8 /xsec. A 50 MHz flash ADC readout system then saves 
the signal traces for all counters in the array with more than 0.3 MIP. Two fits are used 
to reconstruct the properties of the cosmic ray. First, a fit to the times that count ers 



were struck, using the modified Linsley shower-shape function (ITeshima et al. 



198J), is 



made to determine the arrival direction and the core position of the event. Subsequently, 



a lateral distri 



(ITakeda et al. 



jution fit, with the same functional form used by the AGASA experiment 



1998 



20031 ). is employed to find 5(800), the density of shower particles at a 
lateral distance of 800 m from the core. The energy is then estimated by using a look-up 
table in 5(800) and zenith angle determined from an exhaustive Monte Carlo simulation. 



3. Aperture Calculation 

In the ultra-high energy regime, computer-time requirements make it impossible to 
follow every particle when simulating showers. An approximation called thinning is used in 
programs like CORSIKA and AIRES to reduce the computational load by only performing a 
small, statistically representative sample of the air shower simulation. Thinned showers can 
be used for simulation of fluorescence detectors because the fluorescence light comes mostly 
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from near the shower axis where the particle density is extremely high, and the fluctuations 
in the signal are dominated by the Poisson nature of fluorescence photon statistics. But for 
surface detectors, which operate far from the shower core, the number of shower particles 
is low and the thinning approximation fails to represent the intrinsic density fluctuations 
within the shower. To simulate the TA SD accurately, we have developed a procedure 



called "dethinning," where we statistical 



its weighted representative fjStokes et al. 



y reg enerate each group of thinned particles from 



2012|). 



The Monte Carlo simulation of the TA SD has the goal of making an accurate 



representation of the data and our detectors. We start with a 
by the CORSIKA program usin g QGSJET-II-03 ( 



hadronic interactions, FLUKA ( Ferrari et al. 



low-energy hadronic interactions, and EGS4 (jNelson et al. 



Ostapchenko 



2005 



ibrary of showers generated 
to m odel high-energy 



2006 



Battistoni et al. 



20071 ) to model 



19851 ) to model electromagnetic 



interactions. For this libra ry, proton showers were used exclusively because both the 



HiRes composition results ( lAbbasi et al. 



result ( iTameda et al. 



20101 ) and the preliminary TA composition 



20111 ) are consistent with protons generated by QGSJET-II-03. A 
complete representation of calibration and ontime for each surface counter as a function 
of time is also included. Events are then chosen from our sho wer libraries accord ing to 



the spectrum previously measured by the HiRes collaboration ( lAbbasi et al 



2008|). Direct 



comp arisons between data and Monte Carlo show that the result closely resembles the 



data (lAbu-Zayyad et al 



2012al ). Figure [2] shows a comparison of the S(800) of cosmic ray 



showers. The excellent agreement between the data and simulation exemplifies the accuracy 
of our simulation and the resulting efficiency calculation of the SD. 



The selection criteria employed in our analysis are as follows: 



1. Each event must include at least five counters. 

2. The reconstructed primary zenith angle must be less than 45°. 
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3. The reconstructed event core must be more than 1200 m from edge of the array. 

4. Both the timing and lateral distribution fits must have x 2 /degree of freedom value 
less than 4. 

5. The angular uncertainty estimated by the timing fit must be less than 5°. 

6. The fractional uncertainty in 5(800) estimated by the lateral distribution fit must be 
less than 25%. 

Between May 8, 2008 and May 20, 2012, 13,100 events above 10 18 - 2 eV were collected that 
satisfy these criteria. Figure |3] shows the efficiency of reconstruction calculated from the 
TA SD Monte Carlo Program. The values of aperture and exposure for this data set, 
corresponding to the 100% efficiency region, are 920 km 2 sr and 3690 km 2 sryr, respectively. 
For energies above 10 18 2 eV (where the efficiency is ~ 10% of its plateau value) we can 
accurately simulate all air showers, both well- and poorly- reconstructed. The resolution of 
the TA SD energy determination is better than 20% above 10 19 eV. 

The uncertainty in energy scale of the Monte Carlo simulation of an SD is large, 
and possible biases associated with the modeling of hadronic interactions (e.g., from 
extrapolations of cross sections measured at much lower energies) are difficult to determine. 
However, the energy scale uncertainty is experimentally well-controlled for a fluorescence 
detector (FD) since the energy measurement is calorimetric. We therefore correct our energy 
scale to the TA FD using events seen in common between the FD and SD. The observed 
differences between the FD and SD events are well described by a simple proportionality 
relationship, where the SD energy scale is 27% higher than the FD. Figure H] shows a scatter 
plot of FD vs SD energies, where the latter have been rescaled. Events from all three FD 
stations were included in this plot. The two southern FD stations were calibrated using 
independent techniques from the northern station, which consists of reconditioned HiRes 
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fluorescence telescopes. The resulting energy scales are consistent for all TA fluorescence 
detectors. 



4. Spectrum 



Figure [5] shows the spectrum measured by the TA SD, where the differential flux, 
J(E) = d 4 N(E) /dEdAdQdt is multiplied by E 3 , and plotted against logi £. The ankle 
structure and the suppression at the highest energies are clearly visible. A fit to a broken 
power law (BPL) determines the energies of these features. The fit finds the ankle at an 
energy of (4.6 ± 0.3) x 10 18 eV and the suppression at (5.4 ± 0.6) x 10 19 eV. The power 
exponents for the three regions (below the ankle, between the breaks, and above the 
suppression) are —3.34 ± 0.04, —2.67 ± 0.03, and —4.6 ± 0.6 resp ectively. Also shown in 



Figure [5] are t 



re spe ctra reported by AGASA (ITakeda et al. 



( lAbbasi et al 



2003), HiRes (monocu 



ar mode) 



2010bf ). The 



20081 ). and PAO (combined hybrid and SD) (j Abraham et al. 
HiRes and TA SD spectra agree very well, both in the energy region above 10 18 85 eV where 
the TA SD is 100% efficient, and also at lower energies where TA employs a substantial 
efficiency correction. 

A linear extrapolation of the power law below the suppression predicts 58.6 events 
above the break; whereas TA observed only 21 events. This difference corresponds to 
a Poisson probability of 1.44 x 10~ 8 , or 5.5 standard deviations significance. A related 
observable, Ei/ 2 , is the energy at which the integral spectrum falls to 1/2 of its expected 
value in the absence of the GZK cutoff. Under a wide range of assumptions about 
the spectrum of extrag alactic sources, E\/2 is predicted to be 1 19 72 eV for proton s 



teerezinsky et al. 



20061 ). HiRes reported log 10 E = 19.73 ± 0.07 flAbbasi et all 120081 ) . We 



measure log 10 E = 19.72 ± 0.05. 



-11- 



This 5.5 standard dev iation observation p rovides independent confirmation of the GZK 

Furthermore, the energy of the cutoff is 



cutoff observed by HiRes (lAbbasi et al 



2008 



consistent with the interpretation that the composition is protonic. 



Abu-Zayyad et al. 



( 12012al ) includes a description of systematic uncertainties in the SD 
spectrum measurement. The largest source of systematic uncertainty in the spectrum is 
that of the energy scale. Since the SD energy scale is fixed to that of th e TA fluorescence 



detec tors, we take the systematic uncertainty in the SD energy to be 22% ( lAbu-Zayyad et al. 



201ll ). the same as the FD. This propagates into a 37% uncertainty in the flux. We estimate 
the systematic uncertainty in the aperture calculation by removing the event selection 
criteria, one by one, and measuring the ratio of the number of events in the data and in the 
Monte Carlo simulation. This ratio does not change by more than 3% in any energy bin 
above 10 18 ' 2 eV, so we assign this value to be the systematic uncertainty in the aperture. 



5. Conclusions 

We have measured the spectrum of cosmic rays in the energy range 10 182 — 10 20 3 eV 
using the surface detector of the Telescope Array experiment. In the analysis, we have 
introduced a technique, new to the ultra-high energy regime for surface detectors, of 
calculating the surface detector aperture using Monte Carlo simulation, which allows us 
to measure the spectrum even when the SD efficiency is less than 100%. This technique 
includes a dethinning process that enables the simulation of air showers with excellent 
detail. We found that the energy scale of the SD determined from simulations can be 
reconciled with the calorimetric scale of fluorescence detectors by a simple renormalization 
of 27%. 

Two features are seen in the spectrum, the ankle and the high-energy suppression. 
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Fitting the spectrum to a broken power law shows a definite break at an energy of 
(5.4 ± 0.6) x 10 19 eV, which is consistent with the GZK cutoff energy expected for protons. 
An extended spectrum beyond the GZK energy is ruled out with a statistical significance 
of 5.5 standard deviations. Our result is in excellent agreement with that of the HiRes 
experiment where fluorescence detectors were used. This result demonstrates, contrary to 
the AGASA claim, that there is no difference between measurements of the cosmic ray 
spectrum using a fluorescence detector and a surface scintillation array once the energy 
scales are normalized. 

In summary, by combining the two techniques of surface detectors and fluorescence 
detectors (used by the AGASA and HiRes experiments), we have now obtained a consistent 
energy spectrum for ultra-high energy cosmic rays from both techniques. The spectrum 
obtained by our experiment demonstrates spectral features, a dip and a cutoff, consistent 
with the interaction of extra-galactic protons with the cosmic microwave background 
(GZK process). Finally, if we account for a 20% systematic difference in energy scale, our 
measurement is in good agreement with the spectrum reporte d by PAQ with one exc eption: 



2010H ): even 



the GZK break is reported at (2.9 ± 0.2) x 10 19 eV by PAO ([Abraham et all 
with a 20% energy scale correction the difference between the TA and PAO measurements 
is three standard deviations. 
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Fig. 1. — Layout of the Telescope Array experiment. This figure shows the 507 surface 
detector counters deployed on a 1.2 km grid. Fluorescence telescopes overlook the surface 
detector at three sites called Middle Drum (MD) on the north, Long Ridge (LR) on the 
southwest, and Black Rock (BR) on the southeast. The Central Laser Facility (CLF) is 
situated in the center of array, equidistant from all three fluorescence detectors. 
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Fig. 2. — Data and MC comparison of the event 5(800) distributions. The reduced y 2 is 
1.06, indicating a good agreement. 
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Fig. 3. — Efficiency as a function of energy. Both trigger and reconstruction effects are 
included. 
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Fig. 4. — Energy comparison between the TA SD and FD after the 27% normalization has 
been applied to the SD. 
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Fig. 5. — Cosmic ray flux multiplied by E 3 . Solid line shows the fit of the TA SD data to a 
broken power law. 
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